High energy leptons from muons in transit 
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The differential energy distribution for electrons and taus produced from lepton pair production 
from muons in transit through materials is numerically evaluated. We use the differential cross 
section to calculate underground lepton fluxes from an incident atmospheric muon flux, considering 
contributions from both conventional and prompt fluxes. An approximate form for the charged 
current differential neutrino cross section is provided and used to calculate single lepton production 
from atmospheric neutrinos. We compare the fluxes of underground leptons produced from incident 
muons with those produced from incident neutrinos and photons from muon bremsstrahlung. We 
discuss their relevance for underground detectors. 



I. INTRODUCTION 

Atmospheric muon fluxes at sea level have pushed neu- 
trino experiments underground or under ice, however, 
even at a depth of more than 1.5 km under ice, the muon 
to muon neutrino induced muon ratio is on the order of 
10 6 . Much of the work in detectors like AMANDA |l| 
and IceCube [2| has concentrated on upward — > fi 
conversions to avoid the downward muon background. 
The new DeepCore module Q, with a threshold of ~ 10 
GeV and situated within IceCube, is designed to allow 
some of the downward muon neutrino flux to be studied 
where the IceCube detector acts as a muon veto for the 
background to muons produced in the DeepCore detec- 
tor. The DeepCore module also offers the potential to 
measure v T — >■ r and v e — > e conversions [3|. 

The downward atmospheric muon flux itself [3-@] 
presents an opportunity to explore the predictions of 
quantum electrodynamics. A well known process is 
charged lepton pair production [§4T3j. e.g., 

[iA — > /ie + e~ X . 

This process is the primary contributor to the muon en- 
ergy loss, along with bremsstrahlung, photonuclear and 
ionization interactions in transit through materials fl4j . 
Most of the energy loss in electron positron pair pro- 
duction occurs through small changes in muon energy. 
Nevertheless, occasionally the charged lepton produced 
carries a large energy. Given that there are so many 
downward-going muons, one may have the opportunity 
through rare high energy lepton production to explore 
physics at high energy scales. 

Through charged lepton pair production, one has the 
potential to see evidence for the onset of charm parti- 
cle production of muons in the atmosphere [15l - [l7j . At 
high energies, the long decay length of light mesons does 
not allow them to decay before the surface of the earth. 
Studying charged lepton signals underground presents 
another opportunity for determining the contributions 
from charm production to the atmospheric muon flux. 

Charged tau pair production is also of interest. At- 
mospheric production of tau neutrinos is quite low, since 
tau production comes from D s and b meson production 



in cosmic ray interactions with the air nuclei 0, 19]. At 
high energies, neutrino oscillations for — > v T are also 
suppressed. Atmospheric muon production of tau pairs 
may produce a tau signal (accompanied by a muon) or 
ultimately a tau neutrino flux in an energy regime where 
few are expected. 

In Sec. II, we describe how our prior work for calcu- 
lating the cross section for lepton pair production from 
muons scattering with atomic targets can be applied to 
numerical calculations of the energy distribution of the 
leptons produced for fixed muon energy || . The approxi- 
mate analytic form of Tannenbaum [20J can be translated 
to an electron energy distribution which agrees well with 
our numerical results. The tau energy distribution must 
be evaluated numerically. To compare the fluxes of un- 
derground leptons produced from incident atmospheric 
muons with those produced from incident neutrinos, we 
present an approximate form of the neutrino-nucleon dif- 
ferential cross section that agrees well with numerical 
evaluations for a large range of incident neutrino ener- 
gies and energy transfers to the produced charged lep- 
ton. We review in Sec. Ill the steps to go from a muon 
or neutrino flux to an electron or tau flux, and we de- 
scribe the parameterizations for atmospheric muon and 
neutrino fluxes used here. 

Our calculation can be applied to many detector ge- 
ometries. In Sec. IV, we have focused our calculation 
of underground electron production for the large under- 
ground Cherenkov detector IceCube. Because electro- 
magnetic showers produced by electrons are difficult to 
distinguish from those produced by photons, we compare 
the flux high of energy electrons produced via pair pro- 
duction with that of high energy photons produced by 
muon bremsstrahlung. In addition to electron signals, 
we also calculate the fluxes of underground taus produced 
by incident muons and tau neutrinos, and we explore the 
possibility of r pair production signals at IceCube. In 
addition to IceCube, we also explore the possibility of 
taus produced from incident muons and tau neutrinos 
in the mountains surrounding the High Altitude Water 
Cherenkov (HAWC) detector [U. HAWC has the po- 
tential to see the subsequent tau decay showers. In the 
final section, we discuss our results. 
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II. PRODUCTION OF CHARGED LEPTONS 

With our focus on lepton pair production by muons in 
transit, we begin our discussion with our results for the 
electron or tau energy distribution for a fixed incident 
muon energy. Another source of electrons and taus, in 
this case single leptons, is from neutrino charged current 
interactions. We discuss this below. 



A. Lepton pair production 

We have presented in Ref. Q the formulas to evalu- 
ate the differential cross section for a charged muon to 
scatter with a nucleus of charge Z and atomic number 
A. These formulas extend the work of Kel'ner |9| and 
Akhundov et al. [l2[ and others [l3| so that they appli- 
cable to both electron-positron pair production and tau 
pair production. We have not done the full calculation 
required for fi + fi~ production because of the complica- 
tion of identical particles in the final state. 

As we showed in Ref. j8j, for electron positron pair 
production, the differential cross section is dominated 
by very low momentum transfers to the target. Conse- 
quently, coherent scattering with the nucleus dominates 
the cross section, and the static nucleus approximation 
of Kel'ner Q is quite good. In the numerical work below, 
we use the differential cross section as a function of E e , 
the energy of the positron or electron. Our numerical 
results for the electron energy distribution for incident 
muon energies of 10 3 GeV, 10 6 GeV and 10 9 GeV are 
shown by the solid lines in Fig. Q] 

While we use our numerical results for the differential 
electron energy distribution, we note that Tannenbaum 
in Ref. [2(| has an approximate expression for the differ- 
ential cross section as a function of v = (E^ — E'^/E^ 
for an initial muon energy E^ and final muon energy E'^ . 
In the limit of low momentum transfer to the nucleus, 
E^v ~ E e + E s , the sum of energies of the charged lep- 
ton pair. To first approximation, E e ~ vE^/2. In terms 
of v, Tannenbaum has [20] 



da 28 
v— = —Z{Z 



dv 



9?r 

xf(e,v) 



l)(anf (1 



z 2 )ln(l + l)-l 



(1) 



Here 77 is the classical radius of the lepton produced and 
z — vm^/Ami. The function f(e,v) can be expressed in 
two limiting regions as 



ln( vE^/6.67mi 



unscreened 



f(e,v) = { } ~ J x (2) 

lnl 184.15m;/m e Z ) fully screened . 



In comparing our numerical result to those found using 
Eq. ((T|), we have taken the lower of the two values of 
f(e,v). We have modified v to access the high energy 
tail of the distribution, and with v c± 2E e /{E^ — \E e ), 



we can reproduce our numerical results as long as v < 4. 
This gives 



E K = 



2v 



-E u 



v < 4 



(1 - m M /£„) 
(1 + m M /£ M ) 



(3) 



neglecting m e compared to m M . The dashed lines in Fig. 
Q]show the differential cross section for muon production 
of an electron of energy E e for E^ — 10 3 GeV, 10 6 GeV 
and 10 9 GeV with the approximation of Tannenbaum 
using Eq. ©. 




E e [GeV] 

FIG. 1: Differential cross section as a function of electron en- 
ergy for fiA — > jj,e + e~X for fixed muon energy. Here A — 14.3 
and Z = 7.23 for water. The solid line shows our numeri- 
cal result from Ref. [1] and the dashed lines show the ap- 
proximation of Tannenbaum [2(| using v = (E^ — E'^/E^ ~ 
2E e /(E„ - E e /2). 



The approximate form for the differential energy dis- 
tribution of the produced electron does not have an easy 
correspondence to tau pair production [8j. For tau pairs, 
there is a significant contribution from inelastic scatter- 
ing of the muon with the target and higher momentum 
transfers need to be taken into account. For tau pairs, 
we can only use our numerical differential cross section. 
Fig. [2] shows the scattering contributions to the differ- 
ential distribution of the r energy for t + t~ production 
by muons in transit through rock with an initial muon 
energy E^ = 10 6 GeV. 

The contributions in Fig. [5] are the coherent scattering 
with the nucleus (coh), deep inelastic scattering (DIS), 
and scattering with individual nucleons (incoh-N) and 
electrons (incoh-e). The curve showing incoherent scat- 
tering with atomic electrons demonstrates the threshold 
behavior for t + t~ production. 

In what follows, we use the numerical differential cross 
sections for I — e or I — t production by /! — > fj£ + (r X , 
where the incident muons are produced by cosmic ray 
interactions in the atmosphere. 
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FIG. 2: Contributions to the differential cross section as a 
function of tau energy for fiA — s> /j,t + t~X for fixed muon en- 
ergy, En = 10 6 GeV. Here A = 22 and Z = 11 is used for 
standard rock. Indicated are the coherent scattering contri- 
butions (coh), deep inelastic scattering (DIS) and scattering 
with individual nucleons (incoh-N) and electrons (incoh-e). 



B. Neutrino charged current interactions 

Neutrino production of charged leptons is the primary 
focus of IceCube and other underground experiments. To 
compare charged lepton pair production by atmospheric 
muons to neutrino production of a single charged lepton 
in interactions with nucleons, it is helpful to have an 
approximate differential cross section for charged current 
neutrino interactions. 

In principle, the neutrino charged current cross section 
depends on the mass of the charged lepton produced. At 
E v ~ 100 GeV, the v T N charged current cross section, for 
isoscalar nucleon N, is about 80% that for v^N [22l. l23j. 
By E v = 1 TeV, the ratio of the v T N to v^N charged cur- 
rent cross sections is about 0.95. We neglect the lepton 
mass corrections in what follows. 

An approximate form for the charged current differen- 
tial neutrino cross section at low energies can be found, 
e.g., in Ref. 0, 

da v 2vn v G 2 E v 
dy 7T 

dap 2m p GfE v 
dy 



0.2 + 0.05(1 - y) 2 J (4) 
^( .05 + 0.2(1 -yf) (5) 



with inelasticity y = 1 — EgjE v . While reasonable ap- 
proximations to the neutrino nucleon differential cross 
section at E v ~ 1 TeV, Eqs. ((4]) and ([5]) do not ac- 
count for the y dependence and energy dependence that 
come from increased contributions from sea quarks as 
neutrino energies increase. They don't reflect the vio- 
lation of Bjorken scaling as the momentum transfer in- 
creases. 

Given our interest in high energies, we calculated the 
neutrino isoscalar nucleon cross section using the CTEQ6 



parton distribution functions [25[ for E v = 50 — 10 12 GeV 
[26l W\ . The numerical results use a small Bjorken x 
extrapolation at very small x according to a power law 



xq(x,Q 2 ) ~ x~ x Hfl|i2|. We use the following parame- 



terization for the differential cross section, 
da cc 2m v G 2 f E t _ noX i 



dy 



9 J (a(E) + b(E)(l-y) 2 ) 



c(E) 



(6) 



For neutrino scattering, we have two energy regimes split 
by neutrino energy E" = 3.5 x 10 4 GeV, with 



a v = 0.19 - 0.0265 (2.214 -log(£; c /£;)) 2 
b v = 0.036 -0.0344 (1.994 -log(£; c /£;)) 2 
C = 2.3 x 10~ 2 E <E" , 



(7) 



and for higher energies 



p 



0.060 (E C /E) Q - 675 
0.169 (E c /E) - 73 
0.66 x 10 p 

1.453(log(£; c )/log(£)) 6 - 24 



(8) 



E > E v r 



Fig. [3] compares the approximate form of the neutrino- 
nucleon differential cross section with Eq. ([B]) and pa- 
rameters from Eqs. ([7][H]) with numerical results. 

Because the differential cross section is dominated by 
the contribution from sea quarks at high neutrino energy, 
we use the same high energy fit for antineutrino-nucleon 
scattering above E v c — 10 6 GeV. The corresponding pa- 
rameters for lower energy antineutrino-nucleon scattering 
are: 



1 2 x 10 p " 

4 log(E) 4 - 05 



(9) 



a 9 = 4.89 x 10" 

p a = -6.31 x 10" 

b p = 0.177 x 10 p " 

p b = -2.78 x 10" 5 log(S) 5 9 

c & = 4.4 x 10~ 3 £°- 32 E <E V C . 

Fig. [4] shows the approximate form for antincutrinos for 
Ey < 10 5 GeV with the parameters from Eq. ©. 

The parameterizations for the y distribution in neu- 
trino and antineutrino scattering with isoscalar nucleons 
are within about 15% of the result using the CTEQ6 
parton distribution functions [251 ] - 



III. FORMALISM FOR UNDERGROUND 
PRODUCTION OF LEPTONS 

A. Electron production 

The probability for a muon to produce an electron via 
pair production as a function of muon energy (E^) and 
emerging electron energy {E l e ), in a depth interval d£, can 
be written as [23, 1301 



fprodO" 



e,E^,E l e 



didE; 



Na da vmr (E^E l e ) 
A 9 dEl 



(10) 
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FIG. 3: Differential neutrino-nucleon cross section defined by 
Eq. [6] with parameters from Eqs. (J3[8j) . The solid lines repre- 
sent numerical results using the CTEQ6 parton distribution 
functions [25j and the dashed lines are our approximate an- 
alytic formula. Fig. (a) is the differential cross section for 
E v < 3.5 x 10 4 with the fit parameters defined in Eq. J7J- 
The curves represent E v = 50, 100, 500, 1000, 5000, 10 4 GeV 
from bottom to top. Fig. (b) is the differential cross section 
with parameters defined in Eq. ((HJ. The curves represent 
incident neutrino energies E v = 10 5 , 10 6 , 10 s , 10 10 , 10 12 GeV 
from bottom to top. 



where Na — 6.022 x 10 23 is Avogadro's number, A is 
the atomic mass, and p is the density of the material the 
muon is traversing. For electrons produced in a detec- 
tor (contained) which begins at depth D and extends to 
depth D + L max . The event rate is given by 



dN 
~7f 



d£ / dE, 



<1E! Jo Je* " di 

dEl5(Ef - El) 



'■'fi+p.y-^ij. 



di \ dEp-f-p 



N A du^r^E^.El) 



dE, 



0.5 

y 



FIG. 4: Differential antineutrino-nucleon cross section defined 
by Eq. (J6| with parameters from Eq. (j9}. The solid lines 
represent numerical results using the CTEQ6 parton distri- 
bution functions [25j] and the dashed lines are the results of 
our approximate analytic formula. The curves represent in- 
cident antineutrino energies of E& = 100, 10 3 , 10 4 , 10 s GeV 
from bottom to top. 



where we identify the initial electron energy El with the 
electron energy detected E[. (This correspondence can- 
not be made for high energy taus.) In the above equation, 
(j)^ (E^ ,D + £ 7 cos 9) is the differential muon flux at depth 
D + I as a function of energy and zenith angle 0. 

The same formalism can be applied to production of 
electrons via charged current electron neutrino interac- 
tions. For production of single electrons or positrons 
from an incident neutrino or antineutrino the produc- 
tion probability has the same form with the replacement 
da pair {E^E l e )/dE l e ->• da cc (E u ,E l e )/dEi. The incident 
differential muon flux at depth also needs to be replaced 
with an incident differential neutrino flux. 



B. Tau production 

Electrons and positrons must be produced in the un- 
derground detector to be observed, however, very high 
energy taus can persist over long distances, although the 
tau loses energy in transit. A tau with E T = 1 PeV has 
a decay length of jct ~ 50 m. The tau's electromagnetic 
energy loss over that distance is governed by f3 T which 
is a factor of m M /m r smaller than the corresponding en- 
ergy loss parameter for muons. In the high energy limit 
(E l T > 1 TcV) with continuous energy loss, 



dE 



-PtE 



(12) 



The quantity z is the column depth. In the constant /3 T 
limit, the relation between tau initial energy E % T and its 
energy E? after traveling a distance i in a material with 
a constant density p is 



x (E^ ,D + £,t 



(11) 



El = E\ exp(-P T p£) 



(13) 
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since z = pi- We approximate j3 T — 8.5 x 10 7 cm 2 /g for 
standard rock [28|, [29| . 

The tau survival probabilit y, a ccounting for tau energy 
loss and its finite lifetime, is [234301 



P SUIV (Ef,E l T ) =exp 



1 



CT T P T p \ E\ 



(14) 



The short lifetime of the tau means that the survival 
probability goes to zero for low energy taus traversing 
any considerable distance. Because of this, we will focus 
our calculation in the high energy limit. 

We note that at high energies where P S urv — 1, a tau 
track without a tau decay will mimic that of a muon with 
lower energy. The electromagnetic energy loss of the tau 
scales with energy according to Eq. (12), so AE/Az ~ 
—/3 T E T . Since E T is not know a priori, a muon with 
E^ ~ PtEt/P^ will show the same AE/Az. Taus that 
don't decay in a detector will be difficult to distinguish 
from lower energy muons. 

The flux of taus produced from the atmospheric muon 
flux entering rock at sea level and which emerge from rock 
of a total thickness D is, including the survival probabil- 
ity of the tau, 



dN 
dE{ 



D nE' T exp(PT P (D-e)) 

dl \ dE\ 

n J El 



dE„ 



Ei! 



x M+A (.E p ,4cos0)- 



dC 



P suxv (Ef,E^)5(Ef - Et exp(-p P (D - £))) 



D 



d( 



,„ N A da pair (E P ,El) 

0>&u — TP 

A dEl 



x exp 



m T /exp(— P T p{D — £)) 



-CTj3 T p 

x cj) ll+ p.(E ll ,t,cos9) 



El 



1 



(15) 



The delta function in the above equation explicitly en- 
forces the energy loss relation found in Eq. (|T2|) . 

Again, the same formalism for the production of tau 
pairs from an incident muon flux can be applied to the 
production of a single tau (antitau) particle from an in- 
cident tau neutrino (antineutrino) by replacing the inci- 
dent muon flux with an incident tau neutrino flux. One 
also needs to replace the differential pair production cross 
section da pa i r /dE T with the differential charged current 
neutrino cross section do~ cc /dE T , namely Eq. ([5]). 



C. Atmospheric lepton fluxes 

Atmospheric lepton fluxes come from cosmic ray in- 
teractions in the atmosphere, producing mesons which 
decay to leptons 0-S Ha - [l8j . For the differential muon 
flux from pion and kaon decay (the "conventional flux" ) 
at sea level, we use the analytical form from Ref. [31JJ 



which can be written as a function of energy and zenith 
angle, (b^ +tl (E p ,d,9), for depth d = at sea level: 



0.175 (GeVcm 2 srs) 
(i^/GeV) 2 - 72 



(16) 



+ 



1 + £ AI cos6>**/i03 QeV 
0.037 

1 + £ M cos 6»**/810 GeV 



Here 103 GeV and 810 GeV are the pion and kaon crit- 
ical energies, respectively, which separate the high and 
low energy contributions to the atmospheric flux. The 
effective cosine, cos 9**, takes into account the spherical 
geometry of the atmosphere and is given by 31] 



cos 9** = 5(0) cos 0* 
S(9) = 0.986 + 0.014 sec 9. 



(17) 
(18) 



The parameterization for cos 9* can be found in Ap- 
pendix A of Ref. 

An additional contribution to the atmospheric flux 
comes from heavy flavor particle production and decay 
in the atmosphere, the so-called prompt flux. For the 
prompt muon flux, charmed meson production domi- 
nates. There are a number of predictions for the prompt 
muon flux [l5T - fl9j . The vertical prompt flux can be pre- 
dicted from a perturbative QCD calculation which can 
be parametrized by [l5| 



logic (E 3 J 

(E,0,9)) 



-5.37 + 0.0191s + 0.156s 2 
0.0153x 3 (19) 



for x = log 10 (£ , M /GeV). The results for atmospheric 
charm using a dipole model evaluation of the cc cross 
section from Ref. [17j give a lower prompt flux prediction. 
The approximate form for the sum of p + p at sea level 
from Ref. [13 is 



2.33 x 10- 6 (GeVcm 2 srs) 1 
(^/GeV) 2 - 53 (l + £ M cos9**/E 



(20) 

where E = 3.08 x 10 6 GeV. Both of the prompt flux 
formulas are the same for the p + p, + , and 
v e + v e fluxes, since to first approximation, the charmed 
mesons decay to electronic and muonic channels with 
equal branching fractions and the energy distribution of 
the muon and the muon neutrino is about the same. 

The atmospheric muon flux at depth in rock or ice 
depends on electromagnetic en ergy losses of the muon as 
it passes through the material [l4 HH • Using the energy 
loss formula 



A = a + PE 
dz 



(21) 



the muon flux at depth d in a material of constant den- 
sity p can be written, to first approximation assuming 
continuous energy losses (cZ), by 



b^E^d^) ~ 0^^,0,9) eMPpd) 



(22) 
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Here, the muon energy at depth d is related to the 
surface muon energy by 



Ef, 



exp(j3pd)Ef l + (exp(/3pd) - l)a//3 



The exponential factor in Eq. (|22|) comes from 



dE; 

dE„ 



exp{j3pd) . 



(23) 



(24) 



The above formulas are valid for continuous energy losses. 
Fluctuations in energy loss have the effect of increasing 
the down-going muon flux underground 32], however, 
these corrections amount to 5-10% enhancements for a 
depth of 1 km w.e. for muon energies between 100 GeV-1 
TeV [3l|. For our evaluations here, we neglect the correc- 
tions to the underground muon flux due to fluctuations 
in energy loss. 

Fig. [5] compares the contributions to the underground 
muon flux at a depth of d = 1.5 km in ice (A = 14.3). 
We use the two-slope parameterization of Ref. [3l| for 
the energy loss, with 



a = 2.67 x 10~ 3 GeVcm 2 /g 
/3 = 2.4 x 10~ 6 cm 2 /g 



for E^ < 3.53 x 10 4 GeV, and 



a = -6.5 x 10~ 3 GeVcm 2 /g 
(3 = 3.66 x 10~ 6 cm 2 /g 



(25) 



(26) 



at higher energies. We show the contributions from the 
vertical conventional flux given in Eq. (|16|) as well as 
both of the prompt parameterizations given in Eqs. f| 19[) 
(upper) and (|20p (lower prompt curve). The prompt flux 
becomes the dominant contribution for muon energies 
E^ > 10 6 GeV. This is due to the fact that charmed 
mesons decay more quickly than pions and kaons. The 
probability for pion and kaon decays introduces a factor 
of 1/E relative to the probability for charm decays in the 
atmosphere for E < 10 8 GeV. 

For the muon produced tau lepton pairs for a detector 
array like HAWC, we need the atmospheric muon flux at 
an altitude of 4.1 km. The high energy atmospheric muon 
flux at this altitude is approximately the same muon flux 
as at sea level. This is because the majority of the muons 
are produced at an altitude of about 15 km [33| . For the 
energies considered here, at altitudes between 15 km and 
4 km, pion and kaon energy loss through interactions 
with air nuclei are favored over meson decays. 

Finally, to compare the electron and tau pair produc- 
tion rates from muons in transit to the rate for single 
electron and single tau production by electron neutrinos 
and tau neutrinos respectively, we need the atmospheric 
electron and tau neutrino fluxes. At the energies consid- 
ered here, E > 100 GeV, the conventional electron neu- 
trino flux is approximately a factor of 135 smaller than 
the conventional muon flux @. For our calculations here, 




10 100 1000 10* 10° 10° 
EjGeV] 



10 



10° 



10° 



FIG. 5: Contributions to the underground muon flux from at- 
mospheric conventional and prompt fluxes. The solid line rep- 
resents the contribution from the conventional atmospheric 
flux given in Eq. (|16l) . The dashed curve represents the con- 
tribution from the atmospheric prompt flux given by Eq. fj 19p 
while the dot-dashed is for Eq. ()20|) . These contributions are 
for a depth of d = 1.5 km in ice in the vertical direction. 



wc use 



1.30 x 10~ 3 (GeVcm 2 srs) 1 



(^/GeV) 2 - 72 
1 



{27) 



1 + £ M cos (9*7103 GeV 



0.037 



l + cos 0**/81O GeV 



The conventional flux of electron neutrinos has an ap- 
proximate 60:40 ratio of v e : 9 e at E v = 1 TeV [5[, 
a ratio we use here for the full energy range. For the 
prompt atmospheric v e + v e flux, we use Eqs. (|19|) and 
(1201) as two representative fluxes. The prompt neutrino 
to antineutrino ratio is 50:50. 

The tau neutrino flux comes from two sources, oscilla- 
tions from conventional neutrinos, primarily u^, and from 
prompt decays of the D s and b mesons and subsequent 
tau decays [l8L Il9l | . The prompt tau neutrino flux from 
Ref. [lH can be written approximately as 

IxlO-'E 05 (Governors)- 1 

^ + ^(E,0,0) - {E /GeVr ( ' 

1 



1 + {E/l x 10 6 ) - 7 + (E/4 x 10 6 ) 



1.5 



The average height of production of atmospheric leptons 
is at an altitude of ~ 15 km (33J. We are considering the 
tau neutrino downward flux or flux at 45° zenith angle 
at energies above ~ 50 GeV. The oscillation of — > v T 
does not contribute significantly to the tau neutrino flux 
over these distances at the energies of interest, so we do 
not include it in our calculation. 
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IV. RESULTS 
A. Underground electrons in IceCube 

We begin with our results for the rate of electron pro- 
duction in large underground detectors. For underground 
electrons, we have focused our calculation on large under- 
ground Cherenkov detectors such as IceCube. We com- 
pare the flux of electrons produced via pair production 
by muons in transit through the detector with the flux 
produced by electron neutrinos or antineutrinos interact- 
ing via charged current interactions. We then show the 
contributions from muon bremsstrahlung which also pro- 
duce an electromagnetic shower. 

Fig. [6] shows the differential flux of electrons and 
positrons as a function of electron energy produced in 
ice between the vertical depths of 1.5 km < d < 2.5 km. 
For muon induced events, we compare the contribu- 
tions from incident conventional and prompt atmospheric 
muon fluxes, labeled fx — > e in the plot. It is impor- 
tant to note that because the production mechanism is 
fiA — > /ie + e~ X , the total number of high energy events 
comes from the sum of the fj, + p. atmospheric flux. Ac- 
companying the electron is a positron, a muon, and pos- 
sibly evidence of the interaction with the target. For 
comparison, we show the vertical conventional fi e 
rate using the Tannenbaum approximation to da pa i r /dE e 
with the dashed line. 

We also show the contributions to the underground 
differential electron flux from an incident flux of elec- 
tron neutrinos and anti-neutrinos, again in the vertical 
direction. For neutrino induced events, there is only one 
charged lepton produced, so the electron flux come purely 
from the incident neutrino flux while the positron flux 
comes from the incident antineutrino flux. Because de- 
tectors like IceCube have no way of measuring charge, 
we have summed the rates electron and positron events 
together to show the total number of high energy events 
produced from incident neutrinos and antineutrinos to 
better compare with the number of events expected from 
atmospheric muon production. 

Roughly speaking, the conventional atmospheric muon 
flux falls approximately as ~ E~ 4 at high energies and 
the prompt flux falls as ~ E~ 3 . With the steeply falling 
muon flux, the production rate of electrons is dominated 
by the high energy tail of da pa i r /dE e . As can be seen 
in Fig. |51 the vertical underground electron flux is dom- 
inated by the conventional atmospheric muon flux for 
electron energies E e < 10 4 GeV. Around E e ~ 10 5 GeV 
the contributions from the prompt flux start to take over. 
This is in contrast to the crossover point for the muon 
flux, which occurs at an energy an order of magnitude 
higher, at ~ I0 6 GeV. The electron energy distribution, 
with the electron accompanied by a muon, may augment 
efforts to measure the onset of the prompt muon flux. 

For the downward electron neutrino and antineutrino 
fluxes, it is not until E e ~ 10 7 GeV that the contribu- 
tions from incident neutrinos start to be comparable to 




EjGeV] 

FIG. 6: The differential underground electron flux scaled by 
the square of the electron energy for electrons produced in 
ice between the vertical depths 1.5 < d < 2.5 km. The solid 
curves represent the electron flux produced by incident verti- 
cal conventional fluxes of muons and neutrinos given by Eqs. 
(|16[) & (|27|) . The dashed curves labeled prompt represent the 
contribution from an initial prompt flux given by Eq. (|19[) 
while the dot-dash curve is from the prompt flux given in Eq. 
(|20[) . The dashed curve following the conventional fj, — ► e 
curve was calculated using the Tannenbaum approximation 

tO da pair /dE e . 



those from incident muons. The crossover between the 
vertical conventional and prompt electron fluxes occurs 
at E e ~ 10 4 GeV. The electron production rate from 
v e — > e is several orders of magnitude lower than from 
fi — y e, except at the highest energies considered here. 

Although our present analysis only considers atmo- 
spheric muon events in the vertical direction, it can be 
generalized to other zenith angles. There are several 
features to consider with increased zenith angle. First, 
the surface conventional atmospheric muon flux increases 
with angle, for example, by a factor ~ 2 for E^ = fO 6 
when the zenith angle increases from 0° — > 60°. The sur- 
face convential atmospheric electron neutrino flux also 
increases with zenith angle. 

At depth, muon and electron neutrino fluxes are af- 
fected differently. Downward neutrinos experience lit- 
tle attenuation. Even at the highest energies, the neu- 
trino interaction length from interactions with nucleons 
is larger than 10 5 km of ice [2q . Muon energy loss from 
the surface to the detector is an important feature. The 
muon flux relative to the surface value reduces according 
to Eq. (22) in the limit of continuous energy loss. The 
net effect is that even though the surface muon flux at 
non-zero zenith angles is larger than the vertical muon 
flux, the flux for a detector at depth D (and slant depth 
d ~ Dj cos 9) is decreased relative to the vertical muon 
flux at depth D. The underground flux of electrons pro- 
duced in IceCube from an incident conventional atmo- 
spheric muon flux at a zenith angle of 30° is reduced to 
~ 80% of the flux produced from muons in the vertical 
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direction for electron energy E e = 10 GeV. For electron 
energy E e — 10 9 GeV, the electron flux produced by at- 
mospheric muons with a zenith angle of 30° is ~ 90% 
relative to the flux produced by vertical muons. 
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FIG. 7: The differential underground electron and photon 
fluxes scaled by the square of the electron or photon energy 
for particles produced in ice (A = 14.3) between the vertical 
depths 1.5 < d < 2.5 km. Fig. (a) shows the total differential 
flux calculated using Eq. (1 1 1 j) . The solid curves represent the 
electron flux produced by incident vertical conventional and 
prompt (Enberg et al.) fluxes of muons. The dashed curves 
show the conventional and prompt /j, — 5* 7 contribution. Fig. 
(b) represents the differential flux calculating by setting the 



upper bound on v e>1 to v™" x = 0.01 {E\ 
curves are the same as in Fig. (a). 
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At IceCube, the electromagnetic showers produced by 
photons look the same as the electromagnetic showers 
produced by electrons. Muon bremsstrahlung in the 
detector is therefore another source of electromagnetic 
showers. In Fig. Eta), we show with dashed lines the 
flux of photons accompanied by a muon from fiA — > fi'yX 
as a function of photon energy. This flux is evaluated 
using the analytic formula of Ref. [34j for the muon 
bremsstrahlung differential cross section. On the scale of 
the figure, the results are not much changed by including 
the more precise scaling of Ref. [35| . 



Above E ~ 20 GeV, the bremsstrahlung contribution 
dominates the electromagnetic signal. This can be un- 
derstood by the characteristic behavior of the energy dis- 
tribution of photons versus electrons. The electron en- 
ergy distribution from muon pair production falls more 
rapidly with E e than dahrem/dE-y ~ E~ , so when con- 
voluted with the steeply falling atmospheric flux, the 
bremstrahlung contribution dominates. The crossover 
between prompt and conventional photon signals moves 
to higher energies than the electron signals, but it is still 
less than the energy of the crossover of the muon flux 
itself. 

Finally, we remark that the muon produced photon 
flux has a different relation between initial muon energy 
and photon energy than the corresponding muon and 
electron energies. This comes from the E 1 scaling be- 
havior which favors i> 7 = E^/E^ > v e — E e /E^ for fixed 
Efj,. If one could correlate the incident E^ to the out- 
going Ej or E e and consider v e ^ < iQ"? = 0.01, the 
electron flux from conventional muons would dominate 
the the photon flux from conventional muons. The dif- 
ferential flux of photons and electrons with this scaling 
restriction is shown in Fig. [TJb). For such a restriction 
in w e ( 7 ) , the electron flux from conventional muons drops 
by a factor of ~ 10 2 , however, the photon flux from con- 
ventional muons would be reduced by a factor of ~ 10 5 . 
The electron signal dominates the bremsstrahlung signal 
by about a factor of 10 with this restriction. Electro- 
magnetic showers from muons still dominate the electron 
neutrino induced electromagnetic showers. 



B. PeV taus at IceCube 

Muon production of high energy electrons and 
positrons may be a background to searches for taus via 
"lollipop" events [Hj]. At these energies, muons do not 
decay, but taus do, yielding an event with a "muon-like" 
track which has a splash of energy from the tau decay. 
The tau track, before the decay, appears as a muon-like 
track. The high energy electron from \xA — > fie + e~X 
could also leave an energy splash, with a continuing muon 
as part of the event configuration. Furthermore, there is 
the potential to produce tau pairs by fiA — > [it + t~ X. 
Some of these taus will appear at the edge of the detec- 
tor, and some of them will decay in the detector. 

To explore tau production by atmospheric muons and 
atmospheric tau neutrinos and antineutrinos, we have 
calculated the vertical tau flux entering the IceCube de- 
tector at a depth of 1.5 km in ice for tau energies between 
1 — 1000 PeV. Fig. [8] shows the different contributions 
to the tau flux from incident atmospheric lepton fluxes. 
At a depth of 1.5 km, the tau events are dominated by 
the charged current production from the incident prompt 
tau neutrino flux given in Eq. (|28|) . In this energy range, 
the contribution to the underground tau flux from an in- 
cident atmospheric muon comes mainly from the prompt 
flux. 
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To convert to lollipop events, one needs to multiply the 
tau flux by the decay probability. Since the decay length 
at 10 6 GeV is about 50 m, the decay probability over the 
1 km of the detector at that energy is Pdecay — 1, and 
it decreases to Pdecay — 0.2 at E T — 10 8 GeV. The flux 
of decaying taus, where the taus are produced by muons 
or tau neutrinos, is quite small compared to the fluxes 
shown in Fig. [5] 
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FIG. 8: Differential tau flux scaled by the square of the final 
tau energy entering the detector at a depth of 1.5 km in ice. 
The lower solid line corresponds to tau production from a 
vertical incident conventional muon flux given in Eq. (I16|l . 
The dashed and dot-dashed curve represents the tau flux from 
a vertical incident prompt flux given in Eq. ()19p and Eq. (|20p 
respectively. The top solid curve is for the tau flux produced 
with an incident prompt tau neutrino flux given in Eq. (|28|l . 

While we have not done a full scale analysis of the 
number of events with decaying taus in the detectors, we 
can use a characteristic 

4> v+9 = 10~ 8 (cm 2 ssrGeV)" 1 (GeV/£;) 2 (29) 

isotropic neutrino flux and look at the relative normaliza- 
tions. At E v = 10 5 GeV, Eq. gives a flux that is a 
factor of about 30 larger than the prompt v T flux used in 
Fig. [SJ The downward tau flux produced near the edge 
of the detector from this E~ 2 tau neutrino flux would 
remain at least a factor of 30 larger than the prompt 
v r — y t contribution shown. It is clear that the prompt 
tau neutrino flux will be quite difficult to see in IceCubc, 
as will be the high energy tau flux from muons in transit. 

C. Tau production for HAWC 

A different geometry for muon production of taus is to 
use a mountain as the muon conversion volume. For suf- 
ficiently high energies of the produced taus, the taus can 
exit the mountainside and decay in the air. The proposed 
HAWC surface array has the potential to measure the tau 
decay air shower. This detector sits in a mountain saddle 



at an altitude of 4.1 km shielded by mountains on two 
sides. 

For our calculation, we have used a zenith angle of 45° 
in the incident flux and 1 km water equivalent distance of 
rock for the incident muons or tau neutrinos in transit. 
Fig. IH1 shows the energy distribution of taus emerging 
from the rock for both incident muon and prompt tau 
neutrino atmospheric fluxes. We also show the contribu- 
tions from an incident prompt muon atmospheric flux. 

As can be seen in the plot, at energies of E T < 10 5 
GeV the dominant contribution to the emerging tau flux 
comes from the conventional atmospheric muon flux. As 
noted above, a tau energy of E T — 10 5 GeV has a de- 
cay length of ~ 5 m. Even at E T — 10 7 GeV, the decay 
length of 500 m may allow the tau decay to be measured 
by HAWC. At energies above E T ~ 10 6 GeV, the contri- 
bution of taus produced by atmospheric muons provides 
about a 20 — 30% contribution to the total "atmospheric" 
tau flux which is dominated by atmospheric tau neutrino 
conversions. As is the case with taus in IceCube from 
atmospheric sources, the event rates in HAWC would be 
quite small. 
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FIG. 9: Differential tau flux scaled by the square of the fi- 
nal tau energy emerging from 1 km water equivalent of rock. 
We use a zenith angle of 45° for our incident fluxes. The 
lower solid line corresponds to tau production from an inci- 
dent conventional muon flux given in Eq. (|16|) . The dashed 
and dot-dashed curve represents the tau flux from an incident 
prompt flux given in Eq. (|19|) and Eq. (|20|l respectively. The 
top solid curve is for the tau flux produced with an incident 
prompt tau neutrino flux given in Eq. (12811 . 



V. DISCUSSION 

Our evaluation of the pair production cross section 
in Ref. @ already showed that Tannenbaum's param- 
eterization of the cross section Q is a valuable shortcut. 
Here, we have also shown that his differential distribu- 
tion for electron positron pair production, with the iden- 
tification of v = (E„ - Eft/En ~ 2E e /(E„ - E e /2), is 



10 



a reasonably good representation of the numerical eval- 
uation of the electron energy distribution. Using the 
approximate form from Tannenbaum for the differential 
cross section yields errors that are < 30% in the calcula- 
tion of underground electron fluxes for electron energies 
10 GeV < E e < 10 9 GeV. 

Underground electron and photon production by atmo- 
spheric muons may aid in our understanding of the atmo- 
spheric lepton flux itself. The crossover point for prompt 
versus conventional sources in terms of the electron en- 
ergy distribution is at a lower energy than the crossover 
point for muons. Comparing Figs. [5] and [5] shows that 
the crossover happens one order of magnitude lower when 
looking at pair produced electrons as opposed to muons 
for a vertical depth of ~ 1.5 km. Fig. [7]shows that for the 
full bremsstrahlung signal, the crossover point is inter- 
mediate between electrons and the incident muons. The 
additional information gained in studying underground 
electromagnetic signals may aid in the determination of 
the correct charm production model. 

The energy threshold for "lollipop" events from tau 
decays in IceCube is ~ 5 PeV. A 5 PeV tau will produce 
> 200 m track length in the detector 36] . Although po- 
tentially difficult to see, the prompt v T —¥ r flux entering 
the detector receives ~ 20% contribution from high en- 
ergy taus produced by prompt muons in transit through 
ice. The conventional /i — » t flux is suppressed by about 
two orders of magnitude in this energy regime relative to 
prompt v r — > t flux. 

When calculating the underground electron or positron 
fluxes from pair production, there should be rare events 
where one of the leptons comes out with a significant 
fraction of the initial muon energy. In a detector like Ice- 
Cube, it may be difficult to identify the accompanying 
lepton, as well as the muon after scattering. This type of 
event could fake a "lollipop" type signal when looking for 
tau neutrino induced tau events. A comparison of Figs. 
6 and 7 shows that the flux of electrons at an energy of 
a few PeV is about a factor of ~ 100 times the flux of 



atmospheric v T induced taus. Muon bremstrahlung con- 
tributions are even larger, potentially adding to a faked 
signal. 

While muon production of tau pairs in the PeV energy 
range is less than the prompt atmospheric v r contribu- 
tion, at lower energies there is the potential for downward 
secondary v T production from conventional muons via 

fiA -> iit + t~X (30) 

T V T X , 

in, for example, the TeV to PeV energy range. In this 
energy range, the taus decay promptly. 

Our calculations have focused on large Cherenkov de- 
tectors, but the same formalism could be applied for cal- 
culating underground lepton rates at the Indian Neutrino 
Observatory (INO). Due to magnetization, INO will have 
the capability of separating the electron (or r~) signal 
from that of the positron (or r + ). This will allow ob- 
servations at INO to determine the energy distribution 
of the charged partner when looking at pair production 
events, something that could be predicted numerically 
with our evaluation of the differential cross section. 

Electromagnetic interactions, in particular lepton pair 
production by muons in transit through materials, are 
interesting in their own right, not just how they affect 
the energy loss of muons. As neutrino telescopes and air 
shower detectors focus on neutrino induced signals, muon 
signals with high energy electrons or taus may provide 
interesting cross checks to neutrino signals. 
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